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ABSTRACT 

In  this  paper  we  consider  the  routing  of  long-haul  aircraft  from  a 
main  base  to  one  or  more  terminal  bases.  For  these  long-haul  markets,  the 
routing  decision  becomes  critical  because  route  profitability  must  be 
evaluated  for  the  extremely  large  number  of  feasible  routes  covered  by  the 
operation.  In  addition,  the  route  selection  task  is  complicated  by  the 
"pickup-and-del  ivery"  characteristic  of  the  problem.  Therefore,  the 
development  of  an  efficient  procedure  for  selecting  good  candidate  routes 
will  facilitate  the  iterative  flight  scheduling  process  and  may  lead  to 
more  profitable  timetables.  We  define  an  aircraft  routing  problem  that 
captures  the  important  profit-generating  factors  (such  as  intercity  traffic 
estimates,  revenues,  operating  costs  and  aircraft  capacities)  in  the  route 
selection  decision.  We  formulate  this  problem  as  a  mixed  integer  program, 
and  develop  a  Lagrangian-based  solution  procedure  that  exploits  the  special 
structure  of  the  problem.  Computational  results  for  several  test  problems 
indicate  that  the  procedure  is  able  to  select  a  small  number  of  profitable 
candidate  routes,  and  provide  good  bounds  that  confirm  the  near-optlmality 
of  the  generated  solutions. 
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1.  INTRODUCTION  AND  LITERATURE  REVIEW 
The  aircraft  routing  decision  is  one  of  the  most  important  components 
in  the  overall  flight  scheduling  process  of  developing  a  profitable 
operational  timetable  of  flights  for  an  airline.  The  flight  scheduling 
process  consists  of  two  phases:  a  schedule  construction  phase  and  a 
schedule  evaluation  phase.  In  the  schedule  construction  phase,  a  set  of 
aircraft  routes  and  the  frequency  of  service  on  each  route  are  first 
determined  to  maximize  the  profit  generated  from  the  operation,  while 
taking  into  account  the  traffic  estimates  and  revenue  for  every 
origin-destination  pair,  aircraft  characteristics  and  operating  costs,  and 
some  operating  restrictions.  The  construction  phase  is  completed  by 
scheduling  departure  times  and  assigning  aircraft  to  match  the  routing  and 
frequency  decisions.  The  resulting  timetable  is  then  examined  by  operating 
personnel  for  feasibility  and  other  cost  and  performance  considerations  in 
the  schedule  evaluation  phase.  Any  desired  improvements  are  then  fed  back 
into  the  construction  phase,  and  a  revised  set  of  routes  and  the  associated 
frequencies  are  determined.  The  flight  scheduling  process  iterates  between 
these  two  phases  until  a  satisfactory  final  timetable  is  obtained 
(Etschmaier  and  Mathaisel  (1984)). 

Early  research  on  the  aircraft  routing  problem  has  focused  on 
applications  where  the  number  of  alternative  routes  is  relatively  small, 
and  mainly  used  linear  programming  as  the  solution  method  (Dantzig  (1963), 
Kushige  (1963),  and  Miller  (1967)).  These  models  represent  frequencies  as 
continuous  variables;  consequently,  the  solutions  are  very  likely  to 
contain  fractional  frequencies,  which  may  result  in  suboptlmal  integer 
solutions  even  after  applying  some  sophisticated  rounding  procedures.  More 
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recently,  researchers  have  started  using  mixed  integer  programming 
formulations  to  model  and  solve  the  aircraft  routing  problem.    For 
instance,  de  Lamotte  and  Mathaisel  (1983)  used  MPSX-MIP  to  solve  problems 
of  moderate  size  faced  by  a  small  short-haul  airline. 

However,  the  literature  is  still  limited  in  terms  of  solution  methods 
for  larger-scale  aircraft  routing  problems  faced  by  long-haul  carriers. 
For  these  airlines,  the  routing  decision  becomes  complex  because  the  large 
number  of  intermediate  cities  covered  by  their  operations  dramatically 
increases  the  number  of  feasible  routes.  Walker-Powell  (1970)  discusses 
the  Port  Linkage  Problem  for  Qantas  Airline's  Kangaroo  Route,  which  covers 
26  cities  from  Sydney  to  London.  The  cities  are  Indexed  so  that 
connections  are  only  possible  in  the  direction  of  increasing  order.  The 
airline  also  imposed  additional  constraints  based  on  operational 
considerations.  For  instance,  to  simplify  the  crew  scheduling  task, 
several  sets  of  cities,  called  "slip  points",  are  selected,  and  every 
aircraft  is  required  to  stop  at  exactly  one  city  in  each  set;  flights  may 
also  terminate  at  one  of  the  slip  points.  Between  two  slip  point  sets,  at 
most  two  intermediate  stops  are  permitted  for  each  flight.  Given  the 
intercity  demand,  unit  revenue  per  passenger  carried,  fixed  operating  cost 
for  each  origin-destination  pair,  and  aircraft  capacity  restrictions,  the 
Port  Linkage  Problem  seeks  a  profit  maximizing  set  of  routes  from  Sydney  to 
London  (and/or  to  some  slip  points).  Etschnaier  and  Richardson  (1973)  and 
Richardson  (1976)  proposed  a  mixed  integer  programming  formulation  for  this 
problem,  and  developed  a  Benders'  decomposition  algorithm  to  solve  it. 
They  reported  computational  results  for  problems  with  up  to  17  cities,  2 
slip  point  sets,  and  at  most  9  stops. 
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The  multi-stop  aircraft  routing  problem  is  related  to  the  pickup-and- 
delivery  routing  problem  because  we  must  consider  both  embarkation  and 
disembarkation  of  passengers  at  each  intermediate  stop.  The  pickup-and- 
delivery  routing  problem  also  arises  in  the  context  of  school  bus  routing 
(Bodin,  Golden,  Assad  and  Ball  (1983))  and  tractor-trailer  routing  with 
mixed  loads  (Ball,  Golden,  Assad  and  Bodin  (1983)).  However,  this  problem, 
which  is  considered  more  difficult  than  the  pickup-only  (or  delivery-only) 
problem,  is  still  largely  unexplored  by  researchers  (Golden  and  Assad 
(1986)) . 

Since  routing  aircraft  is  the  central  issue  in  the  schedule 
construction  phase,  developing  more  efficient  and  effective  route  selection 
procedures  should  result  in  more  profitable  timetables.  In  this  paper  we 
address  the  multiple  (homogeneous)  aircraft  routing  problem  faced  by  long- 
haul  carriers.  Given  traffic  estimates  and  revenues  for  each  origin- 
destination  pair,  aircraft  operating  costs,  and  aircraft  capacities,  the 
aircraft  routing  problem  seeks  a  set  of  good  candidate  routes  from  a  main 
base  to  one  or  more  terminal  bases  in  order  to  naxinize  total  profit.  We 
assume  that,  for  long-haul  operations,  the  cities  can  be  indexed  so  that 
connections  are  only  possible  from  lower  indexed  nodes  to  higher  indexed 
nodes.  This  restriction  on  the  direction  of  connections  is  considered 
practical  and  even  desirable  because  long-haul  markets  usually  cover  very 
vast  geographic  areas  and  the  distances  between  interaediate  cities  are 
relatively  long.  We  note  here  that  the  aircraft  routing  problem  requires 
Joint  consideration  of  demand  selection  (i.e.,  determining  the  proportion 
of  traffic  demand  to  serve  between  each  origin-destination  pair)  and 
routing  of  aircraft  in  order  to  achieve  the  profit  maximization  objective. 
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In  this  context,  the  demand  selection  decision  arises  because  of  limited 
aircraft  capacity  and  differences  in  the  profitability  of  different  origin- 
destination  pairs. 

The  rest  of  this  paper  is  organized  as  follows:  We  first  present  a 
formal  definition  of  the  long-haul  aircraft  routing  problem,  and  formulate 
it  as  a  mixed  integer  program  In  Section  2.  Rather  than  attempting  to 
solve  the  problem  optimally,  we  focus  on  devising  efficient 
optimization-based  procedures  to  select  good  candidate  routes  and  obtain 
tight  upper  bounds  to  evaluate  the  quality  of  these  results.  We  use 
Lagrangian  relaxation  for  this  purpose.  Section  3  describes  our  Lagrangian 
relaxation  scheme  for  the  mixed  integer  programming  formulation.  We 
develop  an  efficient  procedure  that  exploits  the  special  problem  structure 
to  solve  the  Lagrangian  subproblems  (and  obtain  upper  bounds),  and  discuss 
a  Lagrangian-based  heuristic  which  constructs  good  candidate  routes  (lower 
bounds).  In  Section  4  we  present  computational  results  for  several 
randomly  generated  problems  with  up  to  26  cities  for  the  single-aircraft 
case,  and  23  cities  for  the  4-alrcraft  case.  Finally,  in  Section  5,  we 
suggest  possible  extensions  of  the  model. 

2.  PROBLEM  DEFINITION  AND  FORMULATION 
The  aircraft  routing  problem  that  we  consider  involves  routing  a  fleet 
of  aircraft  of  the  same  type  (i.e.,  aircraft  having  the  same  capacity  and 
operating  costs)  from  a  main  base,  through  several  intermediate  cities,  to 
one  or  more  terminal  bases.  As  indicated  earlier,  the  cities  (including 
the  bases)  are  indexed  so  that  connections  are  possible  from  city  1  to  city 
j  only  if  i  <  j.   The  maximum  (forecasted)  intercity  traffic  demand  (number 
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of  potential  passengers)  for  each  origin-destination  (0-D)  pair  (h,k)  (with 
h  <  k)  is  known,  and  is  denoted  as  d^^j^.  However,  the  entire  traffic  demand 
need  not  be  satisfied.  The  airline  earns  a  revenue  of  rj^j^  for  every 
passenger  it  transports  from  origin  h  to  destination  k.  We  assume  that  the 
total  operating  cost  is  separable  by  route  segments;  the  operating  cost  of 
traveling  on  segment  (i.j).  i.e.,  directly  from  city  i  to  city  j,  is  fixed, 
and  denoted  as  f^i-  As  we  demonstrate  later,  our  model  can  also 
accommodate  variable  operating  costs  that  depend  on  the  traffic  volume  on 
each  leg.  For  simplicity,  we  ignore  these  variable  costs  in  our  model 
development.  The  aircraft  routing  problem  consists  of  (1)  determining  the 
number  of  passengers  to  transport  between  each  0-D  pair  (subject  to  maximum 
traffic  demand  constraints),  (ii)  assigning  the  selected  traffic  to  each 
aircraft  (subject  to  aircraft  capacity  restrictions),  and  (iii)  routing  the 
aircraft.   The  problem  objective  is  to  maximize  total  profit. 

Every  airline  imposes  some  additional  operating  constraints  on  its 
aircraft  routes;  we  do  not  include  all  such  detailed  restrictions  in  our 
model.  Instead,  we  focus  only  on  those  essential  factors,  such  as  traffic 
demand,  revenues,  operating  cost,  and  aircraft  capacity,  that  most 
influence  the  underlying  routing  decision  of  a  long-haul  carrier.  However, 
the  analysis  of  our  basic  aircraft  routing  model  should  provide  useful 
insights  and  information  for  the  actual  route  selection  decision.  First, 
since  the  model  captures  some  of  the  most  relevant  profit-determining 
factors,  it  can  serve  as  a  preliminary  screening  procedure  to  assess  the 
market  profitability  of  the  regions  in  which  the  airline  wishes  to  operate. 
Secondly,  the  routes  selected  by  the  model  represent  potentially  profitable 
routes  that  deserve  further  examination  in  the  iterative  flight  scheduling 
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process.  Finally,  the  profit  indicated  by  this  model  can  serve  as  a 
benchmark  for  cost-effect  analysis  of  other  operational  and  service-level 
considerations . 

The  aircraft  routing  problem  is  defined  over  a  directed  network 
G:(N,A)  consisting  of  n  nodes  in  the  set  N,  and  at  most  n(n-l)/2  directed 
arcs  belonging  to  A.  Nodes  of  the  network  represent  the  various 
intermediate  cities  and  bases,  and  arcs  correspond  to  potential  flight 
segments.  We  assume,  without  loss  of  generality,  that  the  airline 
specifies  a  unique  terminal  base.  Otherwise,  we  connect  the  last  terminal 
base  to  the  remaining  bases  using  dummy  arcs  (that  are  directed  towards  the 
last  base)  with  zero  cost.  The  nodes  are  Indexed  from  1  to  n,  with  indices 
1  and  n  denoting  the  main  base  and  terminal  base,  respectively;  indices  2 
to  n-1  denote  the  intermediate  cities.  The  network  contains  a  directed  arc 
from  city  i  to  city  j  only  if  i  <  j,  and  direct  service  from  i  to  j  is 
possible.  With  this  network  representation,  an  aircraft  route  is  a 
directed  path  from  node  1  to  node  n.  We  will  use  (i,j)  to  denote  the 
directed  arc  or  flight  segment  from  i  to  j ;  (h.k)  denotes  the  0-D  pair  with 
origin  city  h  and  destination  city  k,  and  <h,k>  represents  the  "commodity" 
(which  will  be  explained  later)  between  0-D  pair  (h.k).  The  input 
parameters  for  the  problem  are  as  follows: 

n   number  of  cities  (including  the  bases)  covered  by  the  operation, 

V   number  of  aircraft  (of  the  same  type)  at  the  main  base, 

b   aircraft  capacity, 

d^l^  maximum  (forecasted)  traffic  demand  for  0-D  pair  (h,k),  1  <  h  <  k  <  n, 

r^^   unit  revenue  for  0-D  pair  (h,k),  1  <  h  <  k  <  n. 


f|j  (non-negative)  fixed  routing  cost  per  aircraft  on  flight  segment  (i.j) 
€A. 

To  formulate  the  aircraft  routing  problem  as  a  mixed  integer  program 
we  define  three  sets  of  decision  variables: 

Suj^  number  of  passengers  transported  from  origin  city  h  to  destination  city 
k. 

y^^  number  of  aircraft  traveling  on  flight  segment  (i,j), 

hk 
x^^   number  of  passengers  originating  at  city  h  and  destined  for  city  k  who 

travel  on  flight  segment  (i,j),  l<h<i  <j<k<n. 

The  s-variables  represent  the  demand  selection  decision  that 
determines  which  0-D  pairs  will  be  served,  and  how  much  of  the  traffic 
demand  for  these  0-D  pairs  will  be  satisfied.  The  y-variables  (or, 
aircraft  routing  variables)  form  a  set  of  aircraft  routes  from  the  main 
base  to  the  terminal  base.  Finally,  for  each  0-D  pair  (h,k),  we  treat 
passengers  originating  at  city  h  and  destined  for  city  k  as  a  separate 
commodity  (denoted  as  commodity  <h,k>);  the  x^^  variable  denotes  the  amount 
of  commodity  <h,k>  transported  on  segment  (i,j).  Since  traffic  demand  is 
only  a  forecasted  value  representing  the  expected  number  of  passengers  for 
each  0-D  pair,  the  s  and  x  variables  are  assumed  to  be  continuous.  Using 
these  decision  variables,  we  formulate  the  following  multi-commodity 
flow-based  mixed  integer  programming  formulation  [ARP]  for  the  aircraft 
routing  problem: 


[ARP] 


maximize    Ir,,s,,  -  If  y 


subject  lo: 


Shk    <-    ^ 


k  i-1 

I  xf 5  -  Ixf^  = 

j-iM  j-h 


Ixf5  <  by 


+Sf^,  if  i=h 


ly,,  <-  V 


J-2 


i-1 


for  all  [hM], 


0.      if  h<i<k  for  all  <hM>. 

-s^,  if  i=k 

for  all  (i.j). 


Yy    -  Yy    =0 

xf5  .  s,j  >  0 
y    >  0,  integer 
Range  of  indices:     l<h<i<j<k<a 


(1] 

[2] 

P) 

(4] 
(5) 

for  2  <  i  <  n-1.  (6] 

for  all  <hJo,  [i.j],  [7] 

for  all  <ia;>.  (i.j],  and       [8) 
for  all  (i.j].  [9] 


The  objective  is  to  maximize  total  profit,  which  is  defined  in  (1)  as 
the  total  revenue  (the  first  term)  minus  the  routing  costs  (the  second 
term).  Constraint  (2)  ensures  that  the  amount  of  traffic  served  between 
each  O-D  pair  does  not  exceed  the  corresponding  demand.   Constraints  (3) 
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are  the  commodity  flow  conservation  equations.  Essentially,  they  set  the 
net  flow  of  commodity  <h,k>  leaving  origin  h  and  entering  destination  k 
equal  to  the  number  of  passengers  s^j^  transported  between  0-D  pair  (h,k); 
at  every  other  node  the  inflow  equals  the  outflow  for  this  commodity. 
Constraint  (4)  states  that  the  total  number  of  passengers  transported 
through  any  flight  segment  should  not  exceed  the  total  capacity  of  the 
aircraft  traveling  on  that  segment.  Constraints  (5)  and  (6)  model  the 
aircraft  movements.  Constraint  (5)  specifies  that  at  most  V  aircraft  can 
leave  the  main  base.  Here,  we  implicitly  assume  that  each  aircraft  makes 
at  most  one  trip  during  the  planning  horizon.  Constraint  (6)  specifies 
that  the  number  of  aircraft  entering  each  intermediate  city  must  equal  the 
number  of  aircraft  leaving  that  city.  Finally,  we  impose  the  forcing 
constraint  (7)  which  specifies  that  passengers  cannot  travel  on  a  segment 
(i,j)  unless  we  schedule  aircraft  on  this  segment.  Observe  that  constraint 
(4)  also  imposes  the  same  restriction.  Hence,  constraint  (7)  is  redundant 
in  the  integer  formulation  [ARP] .  However,  they  are  not  redundant  when  we 
relax  the  Integrality  restriction  on  the  y  variables;  thus,  including  them 
in  the  formulation  improves  the  linear  programming  relaxation  bound. 

Notice  that,  without  our  assumptions  about  permissible  flight 
directions  (from  lower  to  higher  indexed  nodes  only),  the  routing  problem 
formulation  requires  additional  constraints  to  prevent  subtours.  A  subtour 
is  a  directed  cycle  that  does  not  include  one  or  both  bases;  thus,  subtours 
are  not  feasible  routes.  However,  when  we  maximize  profit  for  a  routing 
problem  that  is  defined  over  a  general  network,  the  model  may  select 
subtours  that  have  positive  net  revenue.  To  prevent  this  phenomenon,  we 
must  add  explicit  subtour  elimination  constraints;  these  constraints 
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increase  the  formulation  size  and  make  the  problem  more  difficult  to  solve. 
Our  node  indexing  assumption  makes  the  network  acyclic,  thus  obviating  the 
need  for  these  additional  constraints  in  our  aircraft  routing  model. 

The  aircraft  routing  problem  can  also  be  formulated  in  several 
alternative  ways.  For  instance,  Etschmaier  and  Richardson  (1973)  formulate 
the  problem  using  two  sets  of  decision  variables:  integral  routing 
variables  and  continuous  single-commodity  variables  that  represent  the 
number  of  passengers  traveling  on  each  flight  segment.  Although  this 
single-commodity  formulation  is  more  compact,  its  linear  programming 
relaxation  (LP)  bound  is  very  loose.  (The  LP  bound  obtained  using  their 
formulation  for  one  small  test  problem  resulted  in  120%  gap  between  the  LP 
and  Integer  solutions;  the  LP  solution  coincided  with  the  integer  solution 
using  our  formulation.) 

The  size  of  the  mixed  integer  program  for  [ARP]  grows  rapidly  as  the 
number  of  nodes  increases.  For  instance,  for  a  20-city  complete  network 
(with  all  possible  directed  arcs  (i,j),  i  <  j),  the  formulation  contains 
190  Integer  variables,  7505  continuous  variables,  and  9425  constraints. 
The  size  is  more  than  doubled  when  we  consider  a  26-city  problem,  for  which 
the  formulation  will  contain  325  integer  variables,  20800  continuous 
variables,  and  24400  constraints.  Mixed  integer  programs  of  this  size  are 
relatively  difficult  to  solve  optimally  using  branch-and-bound  or  other 
enumeration  methods.  In  particular,  the  imbedded  "plckup-and-delivery" 
characteristic  of  the  problem  complicates  the  task  of  evaluating  the  large 
number  of  feasible  routes  (for  instance,  a  20-city  problem  has  262,144 
feasible  routes;  the  number  of  feasible  routes  increases  to  16.777,216  for 
a  26-city  problem)  .    Therefore,  we  focus  on  developing  an  efficient 
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procedure  to  identify  good  (possibly  suboptimal)  demand  selection  decision 
and  candidate  routes  using  a  Lagrangian  relaxation  scheme  which  we 
introduce  in  the  next  section. 

3.  THE  SOLUTION  METHOD 
A  number  of  successful  applications  of  the  Lagrangian  relaxation 
approach  have  been  reported  in  the  literature.  For  instance,  Lagrangian- 
based  algorithms  have  been  used  effectively  to  solve  the  traveling  salesman 
problem  (Held  and  Karp  (1970)),  the  facility  location  problem  (Cornuejols, 
Fisher  and  Nemhauser  (1977)),  and  the  generalized  assignment  problem 
(Fisher,  Jaikumar  and  Van  Wassenhove  (1980)).  Geoffrion  (1974),  Shapiro 
(1979),  and  Fisher  (1981)  discuss  the  Lagrangian  relaxation  approach  and 
its  numerous  applications.  The  method  consists  of  removing  certain 
"complicating"  constraints,  and  incorporating  them  in  the  objective 
function  using  Lagrangian  multipliers.  The  purpose  of  this  relaxation  is 
to  obtain  a  subproblem  that  can  be  solved  efficiently  because  of  its 
special  structure.  For  a  given  set  of  multipliers,  the  relaxed  problem's 
objective  function  value  serves  as  an  upper  bound  (for  maximization 
objective  functions)  on  the  optimal  profit  of  the  original  problem. 

In  our  Lagrangian  relaxation-based  solution  procedure,  we  dualize  the 

hk 

commodity  flow  conservation  equations  (3)  using  Lagrangian  multipliers  Ui 

for  all  <h,k>,  h  <  i  <  k.  Intuitively,  u^  represents  the  imputed  routing 
cost  to  transport  each  passenger  from  origin  h,  through  city  1,  to 
destination  k.   The  relaxed  problem,  denoted  as  [RP] ,  is  given  below. 
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[RP] 


Zrp  =  maximize    I  r^s,,  -  ^f.^y    -   I  I  c[^jX[>5 

(h.k)  (ij3         ^         tij)<hji> 

subject  to  [2],  (4]  -  [9]. 


Where, 


Thk  =  J'hk  +  <^  -  ^  for  a^l  tWc], 

cf^  =  u|*  -  u]'''  for  all  <lik>.  [i.j]. 


For  a  given  set  of  Lagrangian  Bultipliers ,  the  relaxed  problen  [RP] 
further  decomposes  into  a  demand  selection  subproblem  (containing  only  the 
s  variables)  and  a  capacity  allocation/route  selection  subproblem 
(containing  only  the  x  and  y  variables).  The  selected  Lagrangian 
multiplier  values  determine  the  subproblem  objective  function  coefficients, 
namely,  the  "adjusted"  unit  revenues  f^j^  and  the  imputed  intercity 
transportation  costs  c"^.  Observe  that  if  the  application  context  imposes 
variable  operating  costs  (in  addition  to  the  fixed  costs),  these  costs  can 

hk 

be  added  to  the  imputed  cost  c"^.  As  we  show  In  the  next  section  the 
demand  selection  subproblem  can  be  solved  by  inspection.  On  the  other 
hand,  the  capacity  allocation/route  selection  subproblem  can  be  transformed 
into  a  minimum-cost  flow  problem  for  the  multiple-aircraft  case  or  a 
shortest  path  problem  for  the  single-aircraft  case.  Both  these  problems 
can  be  solved  efficiently  using  specialized  network  flow  algorithms. 
Solving  the  demand  selection  and  capacity  allocation/route  selection 
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subproblems  gives  an  upper  bound  to  [ARP] .   To  improve  the  upper  bound,  we 
iteratively  update  the  Lagrangian  multipliers  using  the  subgradient  search 
method.   We  also  apply  a  Lagrangian-based  heuristic  to  construct  feasible 
solutions  that  provide  lower  bounds  to  [ARP] . 

These  different  components  of  the  Lagrangian  relaxation  approach  are 
described  in  greater  detail  in  the  following  sections:  Section  3.1 
describes  both  the  demand  selection  and  the  capacity  allocation/route 
selection  subproblems,  and  their  solution  methods.  Section  3.2  discusses 
the  subgradient  search  procedure  that  we  use  to  update  the  Lagrangian 
multipliers.  In  Section  3.3,  we  present  a  Lagrangian-based  heuristic  that 
constructs  locally  optimal  solutions  to  [ARP]. 

3.1.  Solving  the  Lagrangian  Subproblems 

The  demand  selection  subproblea,  abbreviated  as  [DSS] ,  consists  of  only 
the  s  variables,  and  is  defined  below: 

[DSS] 

maximize     ^  Thk^hk  ^0) 

subject  to: 

0  <  s^  <  dj^         for  all  (lUc].  (11] 

Observe  that  [DSS]  can  be  solved  by  inspecting  the  sign  of  fj^j^.  That 
is.  if  f^  is  positive  for  0-D  pair  (h,k),  this  0-D  pair  is  "selected",  and 
the  entire  demand  for  this  0-D  pair  is  transported;  otherwise,  no  passenger 


14 

for  this  0-D  pair  is  transported.   Therefore,  the  optimal  solution  to  the 
subproblem  [DSS]  is: 


'hk 


^hk-  if  fhk  >  0 
0.    if  fhk  <  0. 


for  all  0-D  pairs  (h.k). 

The  capacity  allocation/route  selection  subproblem,  denoted  as 
[CARSS],  contains  only  the  x  and  y  variables  and  is  defined  as  follows: 

[CARSS] 


minimize         ^rv  ^  +  X  ^  c|%|^  (12] 


subject  to: 


Sy,j  <  V  [13) 


J-2 
n  1-1 


for  2  <  i  <  n-1. 

[14) 

Yj;^)  <  by 

for  all  (i.j]. 

[15) 

^"5  <-  V, 

for  all  <hJo.  [i.j). 

[16) 

x|^  >  0 

for  all  <hJo.  [i.j).  and 

[17) 

y    >  0.  integer 

for  all  [i.j). 

[18) 

Essentially,  [CARSS]  seeks  to  minimize  the  total  routing  costs  (12)  by 
finding  a  set  of  aircraft  routes  (defined  by  (13)  and  (14)),  and  allocating 
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the  capacity  on  each  segment  of  the  route  among  the  "commodities"  <h,k> 
(restricted  by  (15)  and  (16)). 

Observe  that  since  [CARSS]  does  not  have  commodity  flow  conservation 
requirements,  the  optimal  allocation  of  capacity  among  the  commodities  can 
be  determined  separately  for  each  arc.  Suppose  that  there  are  m(<  V) 
aircraft  traveling  on  arc  (i.j).  The  optimal  capacity  allocation,  and  the 
minimum  cost  corresponding  to  this  aircraft  flow  can  be  determined  by 
solving  the  following  continuous  knapsack  problem  [CKP"j]: 

[CKP^j]: 

6l^  =  mmimize    ^  c^h^)  +  mf„  [19] 

subject  to: 

^xf\  <  mb  (20) 

0  ^  ^'^  ^  <Jhk  for  all  <lUo.  [21] 

This  continuous  knapsack  problem  can  be  solved  efficiently  using  the 
following  greedy  procedure:   Sort  the  commodities  in  increasing  order  of 

Cjj.   Consider  each  commodity  <h,k>  in  sequence  from  this  list.   If  c^^  is 

hk 
negative,  set  x^^   equal  to  the  minimum  of  commodity  <h,k>'s  traffic  demand 

(dj^l^)  and  the  remaining  aircraft  capacity.   Stop  if  the  aircraft  capacity 

is  exhausted  or  all  remaining  commodities  have  non-negative  imputed  routing 

cost  Cjj. 

Notice  that  dj^jj,,  which  denotes  the  minimum  routing  "cost"  when  m 
aircraft  travel  on  arc  (i.j),  is  convex  in  n  since  the  commodities  are 
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selected  in  non-decreasing  order  of  c^^.  (As  the  number  of  aircraft  m 
increases  the  incremental  benefit  of  selecting  additional  passengers 

decreases.)   After  computing  d^ ^^  for  all  arcs  (i,j)  and  m  =  1 V,  we 

can  solve  [CARSS]  by  finding  a  minimum-cost  flow  from  node  1  to  node  n, 
with  the  convex  cost  function  d^^^^  on  each  arc.  (Recall  that,  in  [CARSS], 
the  flow  on  each  arc  represent  the  number  of  aircraft  assigned  to  that 
flight  segment.)  If  the  resulting  minimum  cost  is  negative,  then  the 
directed  paths  in  the  optimal  minimum-cost  flow  solution  define  an  optimal 
set  of  aircraft  routes,  with  the  number  of  aircraft  on  each  route  equal  to 
the  flow  on  the  corresponding  directed  path;  otherwise,  no  route  is 
initiated. 

We  next  demonstrate  how  to  transform  the  minimum  convex  cost  problem 
into  an  ordinary  network  flow  problem  over  a  suitably  modified  network. 
Figure  1  gives  a  small  example  illustrating  this  transformation.  The 
original  network  contains  3  nodes  and  3  directed  arcs.  We  set  up  a 
modified  network  for  the  equivalent  minimum-cost  flow  problem  as  follows: 
To  model  the  convex  cost  function  on  each  original  arc  (l,j).  we  introduce 
V  parallel  directed  arcs  from  node  1  to  node  j  with  costs  equal  to  d^.^  and 
(d|j^  -  d^  j  u_i).  respectively,  for  the  first  and  the  m  arcs,  m  =  1 , 
....  V.  Thus,  the  cost  on  each  parallel  directed  arc  from  node  1  to  node  j 
is  the  marginal  cost  of  assigning  an  additional  aircraft  to  segment  (i,j). 
Each  directed  arc  carries  an  upper  limit  on  flow  of  1  unit.  In  addition  to 
the  parallel  arcs,  the  modifies  network  also  contains  a  dummy  node  D,  and 
two  dummy  directed  arcs  (1,D)  and  (n,D).  These  two  dummy  directed  arcs 
have  costs  equal  to  0,  and  upper  limits  on  flow  equal  to  V.  The  ■Inimum- 
cost  flow  problem  is  solved  from  node  1  to  node  D.   If  the  minimum  cost  is 
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zero,  no  route  is  selected;  otherwise,  the  optimal  arc  flows  define  an 
optimal  set  of  aircraft  routes  for  subproblem  [CARSS] . 


insert  Fig.  1  about  here 


For  the  single-aircraft  case,  since  only  one  cost  coefficient,  d^-j, 
is  necessary  for  each  arc  (i,j),  [CARSS]  can  alternatively  be  solved  by 
finding  a  shortest  path  from  node  1  to  node  n  using  d^jj  as  the  length  for 
arc  (i,j).  Observe  that  even  though  arc  lengths  may  be  negative,  the 
shortest  path  problems  can  be  solved  easily  since  the  network  does  not 
contain  any  directed  cycles  (and,  hence,  any  negative  cycles).  If  the 
length  of  the  shortest  path  is  negative,  the  shortest  path  defines  an 
optimal  aircraft  route;  otherwise,  no  route  is  initiated. 

3.2.  Generating  Upper  Bounds  for  [ARP] 

For  any  given  vector  U  of  Lagranglan  multipliers,  the  objective 
function  value,  Zpp(U),  of  the  Lagranglan  problem  [RP]  provides  an  upper 
bound  to  [ARP].  To  find  the  best  upper  bound  we  must  solve  the  following 
dual  problem: 


[DP]  min  Zpp(U) 

U 


Our  method  attempts  to  find  a  near-optimal  solution  to  [DP]  by  using  a 

hk 
subgradient  search  method  to  update  the  multipliers  u^     and  improve  Zpp(U). 

hk 

Since  the  multipliers  u"  serve  as  imputed  variable  costs  for  transporting 

passengers  from  origin  h,  through  city  i,  to  destination  k,  we  initialize 
.hk 
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'-^hk  /  2.        if  i=h 
"i*^  =  I  (^hi  ^  Sik)  /  b,  if  h  <  i  <  k       for  all  (h.k) 
[  ''hk  /  2.        if  i=k. 

where,  g  is  the  shortest  path  distance  from  node  p  to  node  q  using  the 
fixed  costs  f^g  as  arc  lengths.  Effectively,  the  initialization  scheme 
prorates  the  fixed  routing  costs  to  each  passenger  traveling  from  an  origin 
city  h,  through  an  intermediate  city  i ,  to  a  destination  city  k,  for  h  <  i 
<  k.  At  both  the  origin  city  h  and  the  destination  city  k,  we  use  half  the 
unit  revenue  for  0-D  pair  (h.k),  and  treat  it  as  a  negative  variable  cost 
for  transporting  one  passenger  from  city  h  to  city  k. 

Since  our  relaxation  scheme  dualizes  the  commodity  flow  conservation 

equations  (3)  in  [ARP] ,  the  solutions  s^j    obtained  from  subproblem  [DSS] 

— hk 
and  Xj^j  from  [CARSS]  may  not  conserve  flow  at  each  node.   The  excess 

outflow  of  commodity  <h,k>  at  each  node  i,  h  <  i  £  k,  defines  a  subgradient 

hk 
direction  for  changing  the  Lagrangian  multipliers  Uj^  .    Our  solution 

procedure  iteratively  updates  the  Lagrangian  multipliers  using  a  standard 

subgradient  search  method  (see,  for  example,  Held,  Wolfe  and  Crowder 

(1974),  and  Fisher  (1981)),  solves  the  revised  Lagrangian  subproblems,  and 

computes  the  new  subgradients  until  either  the  number  of  iterations  exceeds 

a  prespecified  limit  or  the  step-size  multiplier  value  falls  below  a  given 

threshold.   We  note  here  that  our  relaxation  scheme  does  not  satisfy  the 

Integrality  Property  (Geoffrion  1974),  and  may,  therefore,  provide  tighter 

bounds  than  the  LP  relaxation  of  [ARP]. 

3.3.  A  Lagrangian-Based  Heuristic  for  [ARP] 

The  Lagrangian  solutions,  although  not  necessarily  feasible  to  [ARP], 
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do  provide  useful  information  for  constructing  feasible  solutions.  Our 
Lagrangian-based  heuristic  uses  the  optimal  solution  to  [CARSS]  to 
construct  an  initial  feasible  solution  for  [ARP].  Then,  we  apply  a  local 
improvement  procedure  that  attempts  to  better  utilize  the  aircraft 
capacity.  Recall  that  the  non-zero  y  solutions  obtained  from  subproblem 
[CARSS]  define  either  a  single  route  or  a  set  of  multiple  routes.  In  the 
following  discussions,  we  will  focus  on  the  case  when  multiple  routes  are 
obtained  from  [CARSS]. 

A  step-by-step  description  of  the  heuristic  follows: 

Step  1:  Separate  the  set  of  routes  into  Vq  single  routes  and  number 
them.  Set  m=l .  Recall  that  the  optimal  solution  to  the  equivalent 
minimum-cost  flow  problem  defines  a  set  of  aircraft  routes,  with  the  number 
of  aircraft  on  each  route  equal  to  the  flow  on  the  corresponding  directed 
path.  The  heuristic  first  obtains  Vq  single  routes  by  tracing  through  each 
directed  path  and  identifying  the  flow  on  the  directed  path. 

Step  2:  Obtain  initial  supply  eunounts  from  the  x  solutions  from 
[CARSS] ■  To  get  a  good  initial  supply  amount  s^j^  for  each  0-D  pair  (h.k) 
that  can  be  served  by  the  set  of  routes  obtained  in  step  1,  we  set  s^ 
equal  to  the  minimum  of  the  outflow  of  <h,k>  from  node  h,  the  Inflow  of 
<h,k>  to  node  k,  the  demand  for  0-D  pair  (h,k),  and  the  aircraft  capacity. 
That  is, 


s°,  =  min  {  I  x,j  .  Ix  ,  .  dj^  .  b  ). 

j-h*l  ■'       j-h  ^ 


This   supply   assignment    is   preferable   to  using  the   s   solution  obtained  from 
subproblem     [DSS],     since     this     subproblem's     "all-or-nothing"     solution 
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characteristic  provides  little  information  on  how  much  of  the  maxlfflum 
traffic  demand  dj^j^  to  actually  transport. 

Step  3:         For  the  m    route,  verify  feasibility  with  respect  to 


aircraft  capacity  constraints .  Let  (H.K)^  denote  the  set  of  0-D  pairs 
(h,k)  that  can  be  served  by  the  m  route.  Although  each  individual  s^j^ 
value,  for  (h.k)  6  (H.Kj^j,  satisfies  the  capacity  constraint,  the  sum  of 
the  s^j^  values  on  some  segments  of  the  route  may  exceed  the  aircraft 
capacity.  Step  3  checks  for  this  possible  situation  and,  if  aircraft 
capacity  is  violated,  reduces  the  supply  amounts  to  those  0-D  pairs  that 
result  in  the  smallest  decrease  in  profit,  until  the  capacity  constraint  is 
met  on  each  segment  of  the  route.  Let  s^j^  denote  the  revised  supply  amount 
for  the  0-D  pair  (h,k) 6  (H,K)^. 

Step  4^:  Compute  possible  improvement  in  profit  by  reallocating 
aircraft  capacity.  For  each  0-D  pair  (h.k)  &  (H,K)|||  that  has  unsatisfied 
demand,  the  heuristic  computes  the  marginal  profit  of  reallocating  the 
aircraft  capacity  by  increasing  the  supply  amount  to  this  0-D  pair,  and 
correspondingly  reducing  the  supply  amount  to  other  0-D  pairs,  if 
necessary.   If  profit  improves  for  none  of  the  0-D  pairs,  go  to  step  6. 

Step  5:  Revise  sj^j^  by  reallocating  the  aircraft  capacity  to  the  most 
profitable  commodity.  The  most  profitable  reallocation  in  step  4  Is 
performed.  After  revising  the  corresponding  supply  amounts,  the  heuristic 
returns  to  step  4  to  check  for  other  possible  improvements. 

Step  6:  Xf  "i^Vq.  Stop.  Otherwise,  check  profitability  of  this  route, 
update  {s^j^} ,  increase  ■  by  1,  and  go  to  step  3.  So  far  we  have  obtained 
an  improved  capacity  allocation.  However,  the  entire  route  may  still 
produce  negative  profit.    In  this  case,  we  simply  discard  this  route, 
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reduce  y^j  by  the  amount  of  aircraft  flow  for  each  (l.j)  on  this  route,  and 

set  Sjjj^  to  0  for  all  (h.k)  6  (H.K)„,.   Finally,  the  initial  set  of  supply 

o  1 
hk' 


amounts  {sPj^>  is  revised  by  setting 


^hk  <—  ""^^  <  s^^  -  s^^.  0  )  for  all  (h.k)  €  (H.K)„. 

The  heuristic  returns  to  step  3  to  examine  the  next  route. 

In  our  implementation,  the  heuristic  is  applied  periodically  (once 
every  10  iterations)  to  improve  the  lower  bound  until  the  step-size 
multiplier  reaches  a  prespecified  threshold  value.  Thereafter,  the 
heuristic  is  applied  after  every  subgradient  iteration. 

4.  COMPUTATIONAL  RESULTS 
Since  actual  data  is  difficult  to  obtain  due  to  its  confidential 
nature,  we  tested  the  solution  procedure  on  the  single-aircraft  12-city 
problem  reported  in  Etschmaier  and  Richardson  (1973),  and  on  several 
randomly-generated  problems.  The  test  problems  reported  here  are  not  only 
larger  (in  terms  of  number  of  cities  and  aircraft)  than  the  previously 
attempted  problems  described  in  the  literature,  but  also  more  closely 
reflect  the  size  of  the  problems  that  arise  in  practice.  The  test  problems 
were  categorized  into  15  problem  classes.  For  each  problem  class,  we 
generated  5  different  instances  using  different  seeds  for  the  random  number 
generator.  Problem  classes  1  to  4  are  single-aircraft  problems  with  17, 
20,  23,  and  26  cities,  respectively.  Problem  classes  5  to  8  are  2-aircraft 
problems  with  17,  20,  23,  and  26  cities,  respectively.  Problem  classes  9 
to  12  are  3-alrcraft  problems  with  17,  20,  23,  and  24  cities,  respectively. 


22 

Finally,  problem  classes  13,  14,  and  15  are  4-alrcraft  problems  with  17, 
20,  and  23  cities,  respectively. 

Given  the  desired  number  of  cities  (including  the  bases),  the  test 
problem  generator  first  randomly  locates  the  required  number  of  nodes  on  a 
100  X  100  grid,  and  numbers  the  nodes  in  increasing  order  of  the  Euclidean 
distance  from  the  origin  (0,0).  All  our  test  problems  are  complete,  i.e., 
the  networks  contain  all  the  possible  directed  arcs.  Demand  for  each  0-D 
pair  was  generated  from  a  uniform  (0,100)  distribution.  Finally,  the  unit 
revenue  for  each  0-D  pair,  and  the  fixed  routing  cost  on  each  arc  were 
generated  as  follows: 

^hk  =  ^'^hk  ""  UNIFORM(0,10)  ,       for  all  (h,k). 
f^j  =  30  *  EU^j  +  UNIFORM(O.IOO) ,  for  all  (i,j). 

where  EU  is  the  Euclidean  distance  between  node  p  and  node  q.  The 
aircraft  capacity,  b,  was  set  equal  to  100  for  all  the  test  problems. 

The  solution  procedure  was  coded  in  FORTRAN  and  implemented  on  an  IBM 
3083  computer.  For  the  single  aircraft  case,  we  solved  the  shortest-path 
subproblem  using  the  algorithm  discussed  in  Lawler  (1976)  for  acyclic 
networks;  for  multiple  aircraft  problems,  we  use  a  modified  version  of  the 
NETFLO  code  (Kennington  and  Helgason  (1980))  to  solve  the  minimum-cost  flow 
subproblem.  Our  implementation  does  not  employ  any  special  data 
structures;  with  improvements  in  data  organization,  the  solution  procedure 
could  perhaps  solve  larger  problems. 

Our  implementation  of  the  subgradient  procedure  initializes  the 
step-size  multiplier  to  2.0,  and  halves  the  multiplier  if  the  Lagrangian 
objective  function  value  does  not  improve  for  20  consecutive  iterations. 
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The  solution  procedure  terminates  when  either  the  number  of  Iterations 
exceeds  1000,  or  the  step-size  multiplier  reduces  to  a  value  below  0.01. 
The  local  improvement  heuristic  was  applied  at  the  first  iteration,  and 
once  every  10  iterations  until  the  step-size  multiplier  reached  a  value  of 
0.1;  thereafter,  the  heuristic  was  applied  at  every  iteration.  For  each 
test  problem  we  recorded  (i)  the  Initial  upper  bound  (obtained  by  solving 
the  Lagrangian  problem  with  the  initial  multipliers),  (11)  the  best  upper 
bound  when  the  procedure  terminates,  and  (iii)  the  best  lower  bound  among 
all  heuristic  solutions.  We  use  the  performance  measure  *GAP  =  (UB  -  LB)  / 
LB  *  100%,  where  UB  and  LB  are  respectively  the  best  upper  and  lower 
bounds,  to  evaluate  the  quality  of  the  solutions. 

Table  1  summarizes  the  performance  of  the  solution  procedure  for  each 
of  the  15  problem  classes.  Notice  that  the  algorithm  was  able  to  reduce 
the  gap  between  the  upper  and  lower  bounds  from  an  Initial  value  of  50.93% 
to  a  final  value  of  1.65%  on  average,  and  the  largest  final  %GAP  was  only 
3.22%  among  the  75  test  problems.  These  statistics  show  the  effectiveness 
of  the  subgradient  method  for  improving  the  upper  bounds.  We  also  observed 
the  following  interesting  characteristic  of  the  heuristic  solutions 
generated  by  the  method.  For  every  problem  instance  that  we  tested, 
subproblem  [CARSS]  generated  only  a  few  alternative  routes  for  step-size 
multiplier  values  below  0.1,  and  the  best  feasible  solution  to  each  problem 
instance  always  resulted  from  one  of  the  routes  In  this  small  set.  Hence, 
the  solution  procedure  seems  to  select  a  small  number  of  good  candidate 
routes  from  the  fairly  large  number  of  possible  routes  in  the  problem. 


insert  Table  1  about  here 
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The  figures  for  CPU  times  in  Table  1  show  that  computational 
requirements  grow  as  the  number  of  nodes  and/or  the  number  of  aircraft 
increases.  However,  the  computational  requirement  increases  at  a  slower 
rate  than  the  problem  size.  For  instance,  the  size  of  a  26-city  complete 
network  problem  is  more  than  twice  the  size  of  a  20-city  complete  network 
problem  in  terms  of  decision  variables  and  constraints  in  the  problem 
formulation.  Yet,  on  average  the  computational  requirement  was  only  1.88 
times  for  the  l-aircraft  case  and  1.75  times  for  the  2-aircraft  case. 
Also,  increasing  the  number  of  aircraft  seems  to  increase  the  CPU  time  only 
marginally.  For  instance,  the  computational  requirement  increases  only 
about  383b  on  average  when  the  number  of  aircraft  increases  from  2  to  4. 

Table  2  records  the  breakdown  of  total  CPU  time  among  the  3  submodules 
of  the  solution  procedure.  For  instance,  solving  the  shortest  path 
subproblem  accounted  for  about  1/5  to  1/4  of  the  total  CPU  time  for  the 
single-aircraft  problems;  solving  the  minimum-cost  flow  subproblem 
accounted  for  about  1/3  of  the  total  CPU  time  for  the  multiple-aircraft 
cases.  Finally,  constructing  feasible  solutions  requires  an  increasing 
proportion  of  the  total  CPU  time  as  the  number  of  nodes  and/or  the  number 
of  aircraft  increases.  We  note  here  that  the  modified  NETFLO  code 
constructs,  in  every  iteration,  a  fresh  network  and  an  initial  feasible 
solution  for  subproblem  [CARSS].  In  retrospect,  since  the  underlying 
network  for  the  minimum-cost  flow  problem  remains  the  same  in  each 
iteration,  some  computational  time  may  be  saved  if  we  had  modified  the 
NETFLO  code  to  retain  the  initial  network.  Furthermore,  since  [CARSS] 
generates  only  limited  number  of  alternative  routes  (for  step-size 
multiplier  values  below  0.1),  we  could  use  the  solution  from  one  iteration 
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to  initialize  the  next,  thus  saving  some  additional  computational  time. 


insert  Table  2  about  here 


We  also  tested  our  solution  procedure  on  the  12-city  problem  listed  in 
Etschmaier  and  Richardson  (1973).  To  ensure  that  our  solution  will  contain 
exactly  one  city  in  the  slip  point  set  (cities  7  and  8),  we  assigned  a 
large  fixed  routing  cost  on  arcs  that  skip  these  nodes  (i.e.,  arcs  (i,j), 
with  1  <  i  <  6  and  9  <  j  <  12),  and  on  arcs  interconnecting  the  nodes  of 
the  set  (i.e.,  arc  (7,8)).  The  large  routing  costs  on  these  arcs  make 
direct  flights  between  the  preceding  (1  to  6)  and  following  cities  (9  to 
12)  unprofitable.  Therefore,  in  order  to  serve  customers  originating  in 
cities  1  to  6  and  destined  for  cities  9  to  12,  aircraft  must  stop  at  either 
city  7  or  city  8.  If  an  aircraft  stops  at  city  7,  the  large  routing  cost 
between  city  7  and  city  8  prevents  a  subsequent  stop  at  city  8.  Therefore, 
our  solution  will  contain  exactly  one  city  in  the  slip  point  set  (cities  7 
and  8).  However,  we  relaxed  Etschmaier  and  Richardson's  original 
restriction  that  aircraft  can  stop  at  most  twice  before  and  after  the  slip 
point  stop.  Table  3  compares  our  solution  with  Etschmaier  and  Richardson's 
schedule. 


insert  Table  3  about  here 


Observe  that  by  relaxing  the  maximum  number  of  stops  restriction,  we 
obtained  a  solution  with  $1078  higher  profit.  This  value  is  essentially  an 
opportunity  cost  for  imposing  the  upper  limit  on  number  of  stops.   Notice 
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also  that  the  optimal  route  for  the  original  problem  can  be  derived  from 
our  solution  by  applying  a  drop  heuristic,  which  deletes  one  of  the  first 
three  cities  that  least  decreases  the  profit.  The  result  of  this  test 
problem  illustrates  the  flexibility  and  usefulness  of  our  model  in 
selecting  good  candidate  routes  which  must  satisfy  some  specific  operating 
constraints.  Similar  procedures  can  be  employed  to  account  for  other 
unique  operational  considerations  in  constructing  aircraft  routes  for  a 
specific  airline. 

5.  SUMMARY 

Flight  scheduling  is  one  of  the  most  important  operational  decisions 
for  an  airline.  The  final  timetable  of  flights  is  the  result  of  several 
iterations  of  the  route  construction  and  route  evaluation  phases.  For 
long-haul  carriers,  the  aircraft  routing  problem  in  the  route  construction 
phase  is  very  complex  because  selecting  a  set  of  good  candidate  routes 
requires  evaluating  a  combinatorially  large  number  of  feasible  routes,  and 
incorporating  the  "pickup-and-delivery"  characteristic  of  the  problem. 
Therefore,  the  development  of  an  efficient  procedure  for  selecting  good 
candidate  routes  will  facilitate  the  iterative  flight  scheduling  process 
and  eventually  lead  to  a  more  profitable  timetable. 

In  this  paper  we  have  discussed  a  version  of  the  aircraft  routing 
problem  that  captures  some  of  the  most  important  profit-determining  factors 
in  the  route  selection  decision  faced  by  a  long-haul  carrier  that  operates 
in  a  base-to-base  environment.  Our  model  makes  the  common  assumption  that 
cities  are  visited  in  increasing  order  of  node  indices.  We  presented  a 
mixed  integer  programming  formulation  for  the  multiple  (homogeneous) 
aircraft  problem,  and  developed  a  Lagranglan  relaxation-based  solution 
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procedure  that  exploits  the  structure  of  the  problem.  Computational 
results  show  that  this  procedure  generates  aircraft  routes  and  traffic 
allocations  that  are  within  1-3%  of  optimality  for  several  relatively 
large,  random  test  problems.  The  proposed  analysis  of  the  aircraft  routing 
problem  serves  as  a  preliminary  screening  procedure  to  assess  the  market 
profitability  of  different  regions,  identify  potentially  profitable  routes, 
and  facilitate  cost-effect  analysis  of  other  operational  and  service  level 
policies . 

In  addition  to  the  basic  problem  setting  described  in  this  paper,  this 
model  can  accommodate  some  operating  constraints  that  may  arise  in 
practice.  For  instance,  in  the  previous  section  we  have  illustrated  a 
method  to  modify  the  cost  matrix  in  order  to  ensure  that  each  route  serves 
exactly  one  city  in  the  slip  point  set.  To  comply  with  other  operating 
constraints  such  as  fixed  number  of  stops  between  slip  point  sets,  we  can 
apply  a  drop/add  heuristic  to  the  best  solution  for  the  basic  model.  We 
note  that  although  the  Lagrangian  upper  bounds  are  still  valid  for  the  more 
restricted  problem,  we  expect  the  %GAP  to  increase  because  imposing 
additional  operating  constraints  will  usually  decrease  the  profit,  and 
hence,  the  lower  bound.  As  we  noted  earlier,  our  model  can  easily 
accommodate  variable  operating  costs  (that  depend  on  traffic  volume)  In 
addition  to  the  fixed  cost  per  aircraft  on  each  flight  segment. 
Application  contexts  where  more  than  one  terminal  base  exits  can  be  solved 
using  a  suitable  network  representation.  Round-trip  operations  can  also  be 
modeled  by  replicating  each  node,  i.e.,  adding  a  mirror  Image  of  the 
original  one-way  network.   The  solution  procedure  can  then  be  applied 
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to  the  expanded  network  with  modified  demand,  revenue,  and  operating  cost 
matrices . 

Several  extensions  and  variants  of  our  basic  model  merit  further 
investigation.  First,  while  our  model  can  accommodate  fixed  as  variable 
operating  costs,  all  our  test  problems  assume  only  a  fixed  cost  for  each 
flight  segment.  Additional  test  results  for  problems  with  variable  costs 
might  be  useful.  Empirical  testing  of  alternative  heuristic  methods  and 
multiplier  adjustment  schemes  is  also  worth  pursuing.  Second,  our  model 
assumes  that  all  aircraft  are  of  the  same  type.  This  assumption  enables  us 
to  use  aggregate  routing  variables  (the  y  variables)  in  formulation  [ARP] , 
instead  of  detailed  binary  variables  to  trace  the  route  for  each  individual 
aircraft.  Modeling  heterogeneous  aircraft  fleet  (with  different  operating 
costs  and  capacities),  on  the  other  hand,  will  require  a  detailed 
formulation  that  destroys  some  of  the  special  structure  in  our  model.  Our 
solution  algorithm  might  possibly  extend  to  this  more  complex  model  as 
well.  We  expect,  however,  that  because  of  the  increased  complexity  the 
method's  computational  requirements  will  increase  and  its  performance 
(measured  in  terms  of  *GAPs)  will  deteriorate  for  problems  with  multiple 
aircraft  types.  Finally,  relaxing  the  node  indexing  assumption  (i.e.,  the 
assumption  that  flights  are  permitted  only  in  the  direction  of  increasing 
node  indices)  might  increase  the  model's  applicability  to  contexts  (such  as 
short-haul  routes)  where  the  cities  are  not  completely  ordered.  Again  this 
extension  increases  the  model  complexity  greatly,  and  may  require  different 
solution  approaches.   Chien  (1987)  explores  some  of  these  assumptions. 
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Captions 


Table  1:  Computational  Results  for  Test  Problems 


Table  2:  Breakdown  of  Total  CPU  Time 


Table  3:  Comparison  of  Solutions  for  the  12-City  Test  Problem 


Fig.  1:  Transforming  [CARSS]  into  a  Minimum-Cost  Flow  Problem 
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Modi  fled  Minimum-Cost  Network 


Table  1 


1  Aircraft 


No.  of 

Nodes 

max 

17 
58  .  76 

20 

23 

26 

Initial 

55.86 

62.38 

64  .  14 

%    GAP 

min 

38.54 

46  .  77 

51  .76 

54.  76 

avg 

41  .82 

51  .  32 

54.87 

56  .  58 

Final 

max 

2  .  17 

1  .89 

1  .21 

1  .28 

%    GAP 

min 

0.94 

0.97 

0  .  98 

1  .  04 

avg 

1  .  66 

1  .23 

1  .09 

1  .  18 

CPU  sees 

max 

74  .87 

142  .  53 

148.23 

200.42 

(Total  ) 

min 

61  .  02 

93.67 

134 . 17 

184.44 

avg 

68  .  09 

102.14 

140.48 

192.49 

2  Aircraft 


No  .  of 

Nodes 

max 

17 
56  .  13 

20 

23 

26 

Initial 

54  .  38 

46.23 

58  .  54 

*  GAP 

min 

43  .  87 

50.22 

40.38 

56.27 

avg 

50.76 

52  .  48 

43.23 

57  .44 

Final 

max 

3  .  22 

2.42 

2  .  58 

2.62 

%    GAP 

min 

1  .32 

1  .  82 

1  .39 

1  .83 

avg 

2.11 

2.  12 

1  .86 

2.22 

CPU  sees 

max 

141 .31 

226.71 

289.01 

363  .  14 

(Total  ) 

min 

130.36 

185.99 

249  .  77 

342 .63 

avg 

136  .  78 

202.64 

269  .  77 

355  .  26 

(continued  on  next  page) 
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Table  1 
( continued ) 


3  Aircraft 


No.  of 

Nodes 

max 

17 

20 

23 

24 

Initial 

48  .  66 

52  .  14 

50  .  36 

52  .71 

*  GAP 

min 

45.  27 

43.03 

47.29 

41.15 

avg 

46  .  35 

48.11 

48.91 

45.  14 

Final 

max 

1  .27 

2.04 

1  .66 

2.63 

%    GAP 

min 

1  .  09 

1  .03 

1  .02 

1  .  03 

avg 

1  .  18 

1  .  48 

1  .38 

1  .76 

CPU  sees* 

max 

174 . 38 

243  .  77 

322.38 

341  .  17 

(Total  ) 

min 

142 . 37 

223 . 86 

287.66 

312.38 

avg 

158 .42 

234  .  46 

304.37 

328.79 

4  Aircraft 


N 

o.  of  Nodes 

max 

17 
58.76 

20 

23 

Initial 

62.14 

60.78 

*  GAP 

min 

48.65 

56  .  71 

53.23 

avg 

52.  17 

58.  29 

56  .44 

Final 

max 

1  .57 

2  .09 

2  .  59 

*  GAP 

min 

1  .  02 

1  .54 

1  .83 

avg 

1  .  39 

1  .82 

2.20 

CPU  sees 

max 

209.38 

298.67 

382  .  77 

(Total  ) 

min 

184  .  23 

267  .  14 

338.13 

avg 

196 . 62 

277.94 

355.72 

Initial  *  GAP  =  (initial  UB  -  best  LB)/(best  LB) 
Final  *  GAP  =  (final  UB  -  best  LB)/(best  LB) 


CPU  times  (excluding  I/O)  in  seconds  on  IBM  3083. 
Each  problem  class  contains  5  different  problem  instances 


Table  2 


1  Aircraft 


No.  of  Nodes 


UB^excl.  Sp2) 

SP 

FEAS^ 


17 


67  .  2% 
20  .  6% 
12.2% 


20 


63  .  8% 
21.4% 
14  .  8% 


23 


61  .  2% 
22  .  1% 
16.  7% 


26 


58  .  0% 
23  .  7% 
18  .  3% 


2  Aircraft 


UB(excl.  NETFLO'*) 

NETFLO 

FEAS 


17 

No  .  0 
20 

f_ 

Nodes 
23 

26 

53.2% 
30  .  4% 
16  .4% 

49.9% 
32  .  5% 
17.6% 

47  .6% 
33  .  1% 
19  .  3% 

43.3% 
35  .  0% 
21.7% 

3  Aircraft 


No.  of  Nodes 


17 


20 


23 


24 


UB(excl 

NETFLO 

FEAS 


NETFLO) 


48  .  2% 
31  .  5% 
20  .  3% 


46  .  5% 
31  .  8% 
21.7% 


44  .  3% 
32  .6% 
23.  1% 


40  .  2% 
34  .  5% 
25  .  3% 


4  Aircraft 


No .  of  Nodes 


UB(excl.  NETFLO) 

NETFLO 

FEAS 


17 


44  .  2% 
32  .  6% 
23  .  2% 


20 


39  .  5% 
34  .  8% 

25  .  7% 


23 


38  .  5% 
35  .  1% 
26  .  4% 


1.  UB :  subgradient  iteration 

2.  SP:  shortest  path  subproblem 

3.  FEAS:  feasible  solution  construction 

4.  NETFLO:  m i n i mum- cos t  flow  subproblem 


Table  3 

Solution  from  Etschmaier  and  Richardson  (1973 
Optimal  Route:  1-3-5-8-9-12 
Profit  :  10312 
Passenger  Distribution: 


3 

5 

8 

9 

12 

1 

60 

11 

0 

0 

0 

3 

79 

0 

0 

0 

5 

41 

49 

0 

8 

10 

31 

9 

59 

Solution  obtained  using  Lagrangian  Relaxation  Algorithm: 
Best  Route:  1-2-3-5-8-9-12 

Profit:  11390  ( UB  =  11451.94.  *GAP  =  0.54*,  CPU*  =  16.83) 
Passenger  Distribution: 


2 

3 

5 
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9 

12 

1 

26 

60 

4 

0 

0 

0 

2 

0 

7 

0 

0 

0 

3 

79 

0 

0 

0 

5 

41 

49 

0 

8 

10 

31 

9 
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CPU  time  in  seconds  on  IBM  3083 
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